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Abstract 

Quantum open systems are described in the Markovian limit by master 
equations in Lindblad form. I argue that common "quantum trajectory" 
techniques corresponding to continuous measurement schemes, which solve 
the master equation by unraveling its evolution into stochastic trajectories 
in Hilbert space, correspond closely to particular sets of decoherent (or con- 
sistent) histories. This is illustrated by a simple model of photon counting. 
An equivalence is shown for these models between standard quantum jumps 
and the orthogonal jumps of Diosi, which have already been shown to cor- 
respond to decoherent histories. This correspondence is compared to simple 
treatments of trajectories based on repeated or continuous measurements. 
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I. INTRODUCTION 



Recently a great deal of work has been done in quantum optics on simulations of contin- 
uously measured systems with dissipation, referred to variously as quantum jumps, relative 
state, and Monte Carlo Wavefunction simulations which are examples of a class of 

techniques known as "quantum trajectories." In these techniques, a deterministic master 
equation for the density matrix of an open system is replaced with a stochastic differential 
equation for a pure quantum state. Averaging the solutions over all realizations of the noise 
reproduces the master equation. Such a stochastic pure-state equation is known as an "un- 
raveling" of the master equation. Unraveling is not unique; in general, there can be many 
stochastic equations which average to the same master equation. A single solution for one 
realization of the noise is called a "quantum trajectory." 

In the original conception, the system was assumed to be monitored by continuous mea- 
surements performed on the environment. The information from the measurements "col- 
lapsed" the system density matrix to a pure state, but the randomness of the measurement 
outcomes made this state unpredictable. The effective description of the continuously mon- 
itored system is thus a stochastic differential equation: an unraveling. Different unravelings 
correspond to different measurement schemes. Averaging over all possible measurement out- 
comes reproduces the master equation, which is the best prediction that could be made in 
the absence of any measurements. 

However, even with no knowledge of the environment and its interaction with the system, 
one can formally unravel any master equation into a stochastic differential equation for pure 
states. In this case, the unraveling is simply used as a means of solving the master equa- 
tion. This provides a useful numerical technique, commonly called "quantum Monte Carlo." 
By averaging over many solutions of the quantum trajectory one can find an approximate 
solution to the master equation. A density operator p on a Hilbert space of dimension N 
requires N 2 — 1 real numbers to represent it; this can be computationally prohibitive for 
large N, while a single state (of size 2N — 2) remains practical, even with the necessity of 
averaging over many runs of the stochastic equation. 

Around the same time that quantum trajectories were introduced, the consistent (or 
decoherent) histories formulation of quantum mechanics was developed by Griffiths, Omnes, 
and Gell-Mann and Hartle J71-0]. In this formalism, one describes a quantum system 
in terms of an exhaustive set of possible histories, which must satisfy a decoherence or 
consistency criterion. Histories which satisfy this criterion do not interfere with each other, 
and may be assigned probabilities which obey the usual classical probability sum rules. 

Both quantum trajectories and consistent histories describe a quantum system in terms 
of alternative possible evolutions; they thus bear a certain resemblance to each other. What 
is more, sets of histories corresponding to possible records of a "classical" measuring device 
will always decohere. Thus, there will be a set of consistent histories which correspond to 
the quantum trajectories of a continuously measured system. 

Exactly such a correspondence has been shown between decoherent histories and quan- 
tum state diffusion (QSD), a particular unraveling of the master equation, by Diosi, Gisin, 
Halliwell and Percival fl3| . QSD trajectories were shown to correspond to a set of approx- 
imately consistent histories for a specific choice of projection operators at closely spaced 
intervals of time. Earlier, Diosi had shown that a particular type of branch-dependent 
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decoherent history first examined by Paz and Zurek [15| corresponded to yet another type 
of trajectory, the orthogonal jump unraveling. Below I will show a similar correspondence 
for quantum jumps, and I conjecture that most useful unravelings correspond to some set 
of decoherent histories in an analogous way. I also demonstrate an interesting correspon- 
dence between standard quantum jumps and the orthogonal jumps of Diosi. I have shown 



the correspondence with decoherent histories for a simple model in an earlier paper |L6 
and there has also been work by Ting Yu |T7| from a rather different perspective on the 
relationship of quantum jumps and decoherent histories. In this paper, I will analyze the 
correspondence between quantum trajectories and decoherent histories in general using the 
concept of it generalized records introduced by Gell-Mann and Hartle [0 , which provides 
a unifying framework for both. Recent work on non-Markovian quantum state diffusion 
explicitly includes the correlation between the trajectory and the state of the environment 
fl9| , |20lj ; in that case, the environment state is precisely this kind of generalized record. 

In section II, I review the formalism of quantum trajectories, giving the quantum jump 
unraveling in some detail but also briefly describing the quantum state diffusion and ortho- 
jump unravelings. I then review the formalism of decoherent (consistent) histories, and 
contrast it with quantum trajectories. I point out that the usual quantum trajectory de- 
scription itself assumes a certain consistency condition. 

In section III, I examine a simple model of a photon counting experiment, in which the 
photodetector is represented by a Markovian environment producing rapid dissipation and 
decoherence in a field mode coupled to a quantum system. This representation reproduces 
the effects of repeated Von Neumann measurements, and in the limit of separated timescales 
between the system and environment allows one to derive a Markovian master equation for 
the system alone. This master equation can be unraveled into quantum jump or ortho- 
jump trajectories as described in sections IA and IB above, and we will see that there is an 
equivalence between these unravelings for this model. 

In section IV, I start with the same model of a photon counting experiment, but this time 
enumerate a set of consistent histories representing different photodetection records. We see 
that the probabilities of these histories are the same as the quantum jump trajectories in 
section III, and that these histories do decohere to a high level of accuracy. 

In section V, I show that the model detector of section III produces generalized records, 
the existence of which guarantees that the consistency conditions will be met. The existence 
of generalized records allows the interpretation of quantum trajectories as consistent histories 
even in a case with no measurement devices; I argue this in the case when the outgoing 
electromagnetic field from the system itself serves as a generalized record (though one with 
many incompatible interpretations). 

Finally, in section VI, the results are summarized, and we see that they both generalize 
the notion of quantum trajectories and provide a useful calculational tool for the formalism 
of decoherent histories. 



II. QUANTUM TRAJECTORIES AND DECOHERENT HISTORIES 
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A. Quantum Trajectories 



The systems of interest in quantum trajectories are described by a Lindblad master 
equation in the Markov approximation |21 |, 

p = -i[H 7 p ] + Y, (l mP Ll - hlL m p - l -pLlL m ^j , (2.1) 

where p is the reduced density operator of the system, H is the system Hamiltonian, and 
the {L m } are a set of Lindblad operators which model the effects of the environment. (Note 
that throughout this paper I use units where h — 1.) One of the notable characteristics 
of master equations of the form (|2.1| ) is that they do not, in general, preserve pure states. 
Suppose that the system is initially in the state 

p(o) = IVoXV'ol- 

Over time it will evolve into a mixture which can be written 



pit) = E \Mt))pdMt)\ = n\Mt))(Mt)\), (2-2) 

where M() is the ensemble mean with probabilities 

Pi > o, J2pz = L 

The decomposition fl2~2] ) is generally not unique; there can be many different sets of states 
{1^)} which give the same density operator p. This ambiguity leads to an ambiguity in 
unraveling the evolution, to which we will return shortly. 

Now let us choose a set of states \ip(t, where £(t) is a random process and 

\il>(t,^(t))) obeys a stochastic evolution equation. By averaging over all solutions of this 
stochastic equation for all realizations of the noise £(t), one finds the solution of Q2.2| ) at 
time t. This stochastic equation will be of the form 

\dif)) = \u)dt+\v)d^{t), (2.3) 

where d£,(t) is a stochastic differential variable representing the random process £(t), and 
can include continuous diffusive noise, discrete jumps, or both. The vectors \u) and \v) 
are functions of the state \ip). (In general, there will be several noise terms with different 
d^{t).) A single solution \if)(t, £(£))) of this equation for a single realization of £(£) follows a 
quantum trajectory in Hilbert space. The complete set of solutions for all £(t) constitutes 
an unraveling of the master equation. 

This is a rather formal notion of quantum trajectories, since it is not clear what phys- 
ical significance the states \ip(t, £(£))) have, if any. Given the existence of many different 
expansions ( |2.2| ) for the same density operator p(t), one must also admit the possibility of 
many different possible unravelings. If the choice of unraveling is not unique, what physical 
reality does it have? Fortunately, it is possible to use the abstract idea of unraveling in 
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practical situations where \ip(t, £(£))) does have a physical interpretation. Indeed, this was 
the original motivation for developing quantum trajectories ||. 

The non-unit arity of the master equation is caused by the loss of information from 
the system to the environment, as their interaction produces entanglement between the 
system and environment degrees of freedom. Within the master equation description this 
entanglement is inaccessible; but if one has experimental access to the environment, it is 
possible to extract information about the system by making selected measurements. If the 
measurements are performed with sufficient precision and frequency, and the initial state is 
known, one can determine the exact state of the system at later times. 

If we include the environment in our description, the joint system-environment state 
remains pure: 

|*> = I>k>® \h), (2.4) 

i 

where |eij) and are states in the system and environment Hilbert spaces, respectively. 
For any state \§>) it is possible to find an expression (|2.4|) such that the |(Xj) are all mutually 
orthogonal, as are the this is the Schmidt decomposition ||22|| . However, in general this 
is not necessary; any decomposition (|2.4|) in which the \bi) (but not necessarily the 
are orthogonal will suffice. In particular, given an observable for the environment, we can 
choose the \b{) to be a basis of eigenstates. Then by measuring the environment, the system 
goes into the state |a») with probability |cj| 2 . 

In this picture the system remains in a pure state, but its evolution is no longer deter- 
ministic: it is influenced by random measurement outcomes. This is often referred to as 
conditional (or relative state) evolution: it is conditioned on the measurement record. If one 
averages over all possible measurement outcomes, a particular decomposition of the density 
operator ( |2.2j ) results. This exactly matches the earlier notion of a quantum trajectory, but 
now the unraveling has a clear-cut physical interpretation: the state \ip(t,£(t))) represents 
our knowledge of the system, conditioned on the random outcomes of a sequence of 
measurements, and is the probability of this outcome. 

Note that not all measurements will satisfy the requirement Q2.2|) ! Indeed, only a fairly 
restricted subclass will do so. We discuss this restriction in HE below, and later in section 
V. 

B. Quantum Jumps 

We can make this clearer with a concrete example. Suppose we consider a quantum 
optical system, such as a small cavity with one mirror partially silvered so that radiation 
can escape. The external electromagnetic field represents the environment in this case, while 
a field mode inside the cavity is the system. We describe the system evolution by a master 
equation 

p = —i[H , p] + ^apa) — (7 /2)a'ap — (7 /2)pa) a, (2.5) 

where Hq represents the Hamiltonian of the system mode, a is the annihilation operator, 
and 7 represents the rate of dissipation from the cavity. We will be considering this simple 
master equation for most of this paper. 
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Equation ( |2.5| ) describes the evolution of the system provided that we know nothing 
about the quanta which escape to the environment, or choose to ignore that information. 
Suppose we place photodetectors outside the cavity in such a way that every escaped photon 
is detected. Each detection gives us some information about the state of the system. In this 
case, with perfect detection, the system evolution becomes 

M = -iH^), (2.6) 
interrupted at random times by sudden quantum jumps 

M-aM, (2-7) 

where H e R is the effective Hamiltonian 

H cS = H - 1(7/2)0*0, (2.8) 

and the jumps occur with a probability given by the norm of the state, as we shall see 
explicitly below Equation ( |2.6| ) has the form of the usual Schrodinger equation, but 

with a non-hermitian Hamiltonian. A state that evolves for a time t without jumping is 
given by 

\m)=e- l ^ t mO))- (2-9) 

The jumps represent detections of emitted photons, which trigger a sudden change in 
our knowledge of the system. The presence of the non-Hermitian terms in the effective 
Hamiltonian represents the effect of not detecting any photons on our knowledge of the 
system; a null measurement thus still affects the system. 

This evolution does not preserve the norm of the state. The actual physical state is 



taken to be \ip) = \ip) / \J (ip\i^) , the renormalized state. An actual physical detector cannot 
determine the time of a photon emission with infinite precision; at best, it will determine 
the time within a short interval At. The probability that an initial state evolves for a 
time T and undergoes N jumps at times t\, . . . , tjv (which are assumed to be widely spaced 
with respect to At) is 

p(|<0}) = ^A^^Trje-^ff^-^^e-^ff^^-^a- ■ - ae-^ff' 1 

x |^)(^|e^W t ---aV^ T - t ^j. (2.10) 

This trace is the norm of the state 

\^t u ...,t N ) = e-^^-^ae-^^-^-^a ■ ■ ■ ae' iAcBh (2.11) 

which is the unrenormalized state resulting from this particular sequence of jumps ti, . . . , t N . 
The norm of the unrenormalized state thus gives the probability density for that state to 
be realized. The density operator p is given by averaging over all realizations with 
probability measure (|2.10|), 
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P (t) = Mmt))m)\) 

= £l$p(l$)$l 

= Yl (7<lK^)(4,...,tJ. 
IV'ii,....t JV > 

= £7^ / T rfti--- / T ^l^,..,^)^,..,^!- (2-12) 



The density operator defined by (|2.12 ) solves the master equation 



It is possible to rewrite these equations in explicitly norm-preserving form, 

|#) = -iH eS \ip)dt + { 1 /2)(a)a)\i))dt+ ( = - 1 ] |^>c?iV, (2.13) 



at the cost of a little extra complexity and nonlinearity ||. This is a stochastic differential 
equation, where dN is a stochastic variable which is except at random times (corresponding 
to the jumps) when it becomes 1. It has statistics 

dNdN = dN, M w (dN) =7<a + a), (2.14) 

where Mi^\ denotes the ensemble average over all trajectories which are in state \ifi) at the 
given time t. The nonlinear form has useful properties, but for our purposes the linear form 
will usually prove more convenient. 



C. Other Unravelings 

The quantum jump equation is not the only way of unraveling the master equation into 
quantum trajectories. In fact, there are an infinite number of such unravelings. Certain 
choices are more commonly used than others, however, so we will go over them quickly. 

One of the most useful is the quantum state diffusion (QSD) equation [J23|] 

|^) = -iH\f)dt + £ ((Ll)L m - hlL m - h(L m )\ 2 ) \f)dt + £ {p m - (L m )) |V>Rm- 

(2.15) 

This is an Ito stochastic differential equation |2~4]| , with the d£ m representing continuous 
complex stochastic processes with ensemble means 

M{dC m ) = M(d£ m d£n) = 0, M(d£ m d£) = 6 mn dt. (2.16) 



It is not difficult to show that this equation also obeys the master equation fl2.1|) in the 
mean: 

p = M(|^>(^|). (2.17) 
Rather than discrete jumps, the solutions of this equation undergo continuous diffusion. 
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This standard form of the QSD equation is nonlinear and norm-preserving. However, 
there is a different unraveling known as linear QSD with equation 



|#> = -iH\^)dt - ^Y, L l L ™\^) dt + Y,Lrn\?P)dU- (2.18) 

m m 

This does not preserve the norm of but is also a valid unraveling of the master equation, 
and has properties similar to the nonlinear equation p5*| . 

The QSD equation was discovered by Gisin and Percival, following from earlier work by 
Gisin in the theory of measurement. Carmichael |J showed that a similar diffusion equation 
arose in the case (|2.5| ) from a relative state approach, just as in quantum jumps, but with 
direct photodetection replaced by balanced homodyne detection. It was shown by Wiseman 
and Milburn [[2(J and others [^7| that the exact QSD equation arises in the case of balanced 



heterodyne detection. Thus, one can consider both approaches as giving the state of a system 
conditioned on a measurement record, but using different measurement schemes. We will 
discuss this further below. (The linear equation (|2.18| ) does not have such a straightforward 



interpretation in terms of measurements, though that does not necessarily imply that such 
an interpretation might not be found.) 

Conversely, Gisin and Percival have shown |23| that jump-like behavior can be exhibited 
by the QSD equation by explicitly including a portion of the photodetector in addition to 
the system; their explicit inclusion of part of the environment is similar to the approach of 
this paper. 

Related to QSD are the orthogonal jumps of Diosi pgfl . These trajectories obey an 
equation 

|#> = -iH\if>)dt + £ [(Ll)L m - hlL m + \{Llt m ) - |(L m )| 2 ) \^)dt 

+ J2(\^m)-m)dN m , (2.19) 

m 

where the dN m now represent a stochastic jump process, 

alN m dN n = 5 mn dN m , M(dN m ) = r m dt (2.20) 

and the set of states \ip m ) are mutually orthogonal and orthogonal to The formula for 
these orthogonal states and their jump rates r m is complicated; fortunately, for the purposes 
of this paper we will not need to worry about it. The deterministic portion of equation 
( j2.19j ) is identical to that of the QSD equation expressed in Stratonovich rather than Ito 
form |21|. 

These orthogonal jumps are in some ways the most economical unraveling of the master 
equation, in that these jumps occur fairly infrequently but cause a large change in the state. 
This unraveling also has close ties to decoherent histories, as we shall see. However, relatively 



little work has been done on the connection between ortho-jumps and measurements ||30| . |3l 
An important result is that ortho-jumps result from choosing the Schmidt decomposition of 
( p.4|) , where the |aj) as well as the \bi) are all orthogonal, and measuring the eigenbasis \bi). 



This measurement gives the maximum average information about the system, but requires 
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a dynamic measurement scheme, since the basis \bi) will in general be different at different 
times. 

This form of ortho-jumps is also explicitly norm-preserving, which results in a nonlinear 
stochastic equation, just as for QSD. Unlike the other equations presented, there is no com- 
pletely linear version of this unraveling; the orthogonalization procedure used in determining 
the \ip m ) is intrinsically nonlinear. However, for certain special problems these \ipm) are lin- 
ear functions of In this special case it is possible to find a linear version of ortho-jumps 
which does not preserve the norm, and which is closely related to the standard quantum 
jump equation. We will examine such a case below. 

All of these unravelings correspond to the same Markovian master equation. If considered 
simply as a way of solving for the density operator, it doesn't matter which unraveling one 
picks; averaging over any set of trajectories will give the same final result. But if one actually 
has a physical description of the environment, and of the measurements being performed on 
it, this equivalence is broken, and a particular choice of unraveling is singled out. This is 
the case in this paper, where we assume a particular (albeit simplified) description of the 
environment, including the effects of a measurement apparatus, to model photodetection. 



D. Decoherent (or Consistent) Histories 

Let us now turn to consistent histories. In nonrelativistic quantum mechanics one can 
specify a set of possible histories by choosing a sequence of times £1,^2, • • • ,t n and at each 
time ti a complete set of orthogonal projection operators {V^,.(tj)}, such that 

K^Kdh) = Kfi^Sa^, = '• ( 2 - 21 ) 

where the P's are Heisenberg operators. These projectors represent a complete set of exclu- 
sive alternatives at the given times. A single history corresponds to a choice of one projection 
operator V % a .{ti) at each time t^. This history can be represented by the sequence of indices 
«i, . . . , a n , which I will denote by h. 

This is not the most general type of history. For instance, one can make the set of 
projections at a later time depend on the choices {a{\ at earlier times, making the histories 
branch- dependent. We will not need this extra generality in this paper, but it can be useful 
in describing measurements which are conditioned on the results of earlier measurements. 

The decoherence (or consistency) criterion is described by the decoherence functional 
D[h,h'], a complex functional on pairs of histories. Two histories h and hi are said to 
decohere if they satisfy the relation 

D[h,ti] = p(h)S hh >, (2.22) 

where p(h) defines to be the probability of history h. A set of histories {h} is said to be 
complete and consistent if all pairs of histories satisfy fl2.22| ) and their probabilities sum to 
unity. 

We define a history operator Ch 

Ch = KS^K^n-x) ■ ■■V 1 ai (t 1 ). (2.23) 
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In terms of this operator the decoherence functional becomes 

D[h,ti]=Tr{C hPo 6l,}, (2.24) 
where po is the initial state. If po is pure, po = then D[h, h'} is an inner product 

D[h,h'] = (^\Cl,C h m, (2.25) 

and the consistency condition ( 2.22|) amounts to the assertion that the states {Ch\ip)} are 
orthogonal for different h, with the probability of a history given by the norm of the corre- 
sponding state. 



E. Consistency of quantum trajectories 

It is often convenient to consider these projections in the Schrodinger rather than the 
Heisenberg picture. In this case, the V's are no longer time-dependent, and instead the 
time-evolution is given explicitly: 

(J h — *p n Q—iH{t n —t n -i)'pn—l g— iH(t n ~i— t n -2) . . .-pi g — iH{t\— to) ^ 26) 

Written in this form, it is clear that the state Ch\ip) has a strong resemblance to our previous 
definition of a quantum trajectory. This succession of Hamiltonian and projection terms 
resembles some kind of pure state evolution, and if no single history has probability 1 this 
evolution will have a stochastic component. 

However, there are a number of important differences. The first is the most obvious: in 
our definition of quantum trajectories, there was no requirement that they should obey the 
decoherence criterion. This is related to the fact that the decomposition ( |2.2|) need not be 
in terms of orthogonal states. More on this below. 

The second difference is that quantum trajectories were framed in terms of a split be- 
tween system and environment degrees of freedom, with the environment traced out of the 
equations of motion in the Markovian approximation. These assumptions are not made in 
consistent histories, which can be quite general. 



The open systems case is important, however |32|-|35|, and is worth treating in depth 
in consistent histories. Most such treatments to date |IT , P5 ,S7[ involve choosing a set of 
projection operators on the system alone, while making no assertions about the environment 
degrees of freedom. That is, one specializes to projections of the form 

V = V sys ® lenv, (2.27) 

where V sys is a projection operator in the Hilbert space of the system and l env is the identity 
operator in the Hilbert space of the environment. The Markovian approximation will only 
be valid for certain choices of environments and of the system/environment interaction, 
of course; but if one is treating the same systems in both the quantum trajectories and 
decoherent histories approach, the approximation should be good in both cases. One ends 
up with a reduced description in terms of the system alone, quite in the spirit of the master 
equation. 
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If we physically interpret the quantum trajectories as evolution conditioned on measure- 
ments of the environment, we realize that (|2.27 ) is not adequate to describe this situation. 



To make statements about the state of the environment, we require projectors onto the 
environment, not just the system. In fact, the correct projectors in this case have the form 

^ = lsys®Penv, (2.28) 

quite different from the usual discussion of decoherence in open systems. Thus, to look for 
a set of histories equivalent to a particular unraveling, we must retain in our description 
enough of the environment to include suitable projections of the form ( |2.28| ). As we shall 
see, in some cases this may be very little. To reduce to a description of the system alone, a 
further tracing out of the environment is then required. 

As far as consistency is concerned, it has not been widely appreciated that quantum 
trajectories (as they are commonly defined) must obey a consistency condition of their 
own. This is easily seen by considering the effects of environmental measurements. Rather 
than tracing out the environment to get a master equation ( |2.1| ), we retain a complete 



description of the system and environment (|2.4| ). There is a Hamiltonian H for the system 



and environment together, such that they tend to become entangled with time: 

\V(t)) =e~ i6t \y{0)). (2.29) 

Suppose we now perform a series of n measurements on the environment spaced At = 
t/n apart in time, with outcomes corresponding to projections of the form ( |2.28|) . The 
(unnormalized) joint state conditioned on these outcomes then becomes 

l*C> = ^e-™^ • ■ A|¥(0)>, (2.30) 



with probability = (^|^). We now want to impose the condition (|2.2j ). In order for 
this to hold we must have 



E Tr ««r{l« r €><« r €l} = E Tr -^{l*€>(*cl}» ( 2 - 31 ) 



which implies 



£Tr env {|^>(^|} = 0. (2.32) 



This equation is by no means automatically satisfied. It is, in fact, a consistency condition 
rather like fl2.22|) . In some ways it is much weaker, since ( |2.32| ) is only imposed on sums 



of off-diagonal terms rather than individual terms (although it must be satisfied for all 
times t). In another way it is stronger: the trace in ( |2.32|) is only over the environment, 
so that an entire operator on the system space is required to vanish, not just the trace. 



This is rather like the partial trace decoherence of Finkelstein [38], discussed in section IV. 
In physical terms, (|2.32|) means that the measurements made must not alter the state of 
the environment such that it then acts back and alters the dynamics of the system. As 
we shall see, quantum jumps (like most of the commonly considered unravelings) satisfies 
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both ( |2.32j ) and ([2.22 ) , and in fact satisfies a consistency criterion stronger than either. 
The Born-Markov approximation which is generally made in deriving the master equation 
explicitly assumes that the state of the environment is not changed by interacting with the 
system; so long as measurements are made only on the outgoing field, (|2.32|) will be satisfied. 

What makes a set of histories decohere? While a general characterization is still the 
subject of research, common mechanisms of decoherence have been widely studied in simple 
models, and are now fairly well understood. Suppose that the state of the system at the 
time £j becomes correlated with some degrees of freedom of the environment described by a 
Hilbert subspace Hi, and the projection operator V l a . singles out a particular state of Hi. If 
the state of these degrees of freedom is subsequently unaltered by further interaction with 
the system, then different choices of the V ai at time ti will render the different histories 
orthogonal, and hence consistent, at all later times. 

This gives us the following picture: the total Hilbert space for the system and environ- 
ment has the form H = H sys ® Hi ® • • -Hn, and the system only interacts with a single 
subsystem Hi at a given time U. We choose a set of projections {V^.} onto the state of Hi at 
each time. Clearly, any history composed of such projections must decohere. These subsys- 
tems Hi are called "generalized records" by Gell-Mann and Hartle, because they record the 
choice V l a for all later times; they are "generalized" because the particular set of degrees of 
freedom represented by Hi might be very complicated, and inaccessible to an experimenter 
for all practical purposes NTS . 



This is highly idealized, of course. In a real system it is usually impossible to identify 
these subsystems Hi] the choices of projection are generally not onto these subsystems 
directly, but instead depend on correlations in the state of the system and environment; 
the interaction with the system does not have this precise, time-limited form; nor does the 
system generally interact directly with the generalized record. However, in some cases one 
can explicitly identify physical systems which serve as generalized records (for instance, 
photons which escape to infinity without interacting further); and even where one cannot, 



it is not unreasonable to assume that something corresponding to them is present ||39|| . It 
may be that in a real environment no subsystem remains forever isolated, but it might be 
effectively true on timescales very long compared to the duration of the experiment, or even 
the lifetime of the universe. In particular, situations which correspond to measurements can 
almost always be assumed to have some set of generalized records. Sometimes these are 
easy to identify: lines drawn on graph paper by a measurement apparatus are more than 
sufficient to cause decoherence of the measurement results at all later times. 



III. THE MODEL 

A. Model of photodetection 

Photomultiplier tubes are very complicated devices, with the complication doing little 
to improve understanding of the detection process. Therefore I will use a detector model 
which is simpler to analyze, based on the techniques used by Haroche et al. |4(| in studying 
cavity QED in microwave cavities. Models of this type have been used before by a number 
of authors |}4T| , f42"|| . Figure 1 illustrates the setup. 
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The heart of our detector is a lossy cavity which is probed by a beam of atoms. We 
assume that the atoms are fairly well localized, so that their motion can be treated as 
essentially classical. One mode of this cavity can be excited by photons leaking in from the 
outside; this mode has creation and annihilation operators and b. We assume that the 
electronic states of the atoms are prepared in a superposition of two levels, which we label 
|0) and |1). If there is an atom in state |1) in the cavity it shifts the resonance frequency. 
Cavity loss results from interaction between the mode and the internal degrees of freedom 
of the cavity, which we model as a reservoir of harmonic oscillators. The time-dependent 
Hamiltonian for this detector model is 

H det (t) = wjttb + £ w(t - tjpbllj}^ + $>*(& + &)(r k + f{) + oj k r\r k , (3.1) 
j k 

where Wdet is the frequency of the cavity mode, is a projector onto state |1) of the jth 

atom, f\ and f k are the creation and annihilation operators for the /cth reservoir oscillator, 
uj k is its frequency, and 7^ is the coupling to the cavity mode. The function w(t) is a window 
function peaked around 0, which represents the strength of the interaction when the atom 
is inside the cavity. It is assumed that the traversal times tj of the atoms are spaced an 
average time r apart, and that the width of the window function w(t) is narrow compared 
to t. The important quantity is the integral 



oc 



2A = / w(t) dt. (3.2) 



The effects of the reservoir coupling have been widely studied and are very well under- 
stood. If we assume that the temperature of the reservoir is low compared to the oscillator 
energy c^dot and the frequencies u k are distributed appropriately, the reservoir produces an 
effective dissipation rate I\ |32| , |35 



The incoming atoms are prepared in the state j^in) = (|0) + |l))/\/2- If the cavity mode 
is in the number state \n) the atoms leave the cavity in the state (|0) + exp(— 2zAn)|l))/v / 2. 
This entanglement between the atoms and the cavity mode suppresses interference terms 
between number states n and n' by a factor of cos(A(n — n')) per atom, where we assume 
|A(n— n')\ < 7r/2. With the atoms spaced approximately r apart on average, this corresponds 
to a decoherence rate of T 2 = — (1/r) lncos(A(n — n')). 

We assume that the dissipation time 1/Ti is long compared to the decoherence time 
l/r 2 , so that every photon that enters the cavity persists long enough to be detected, but 
is short compared to the rate at which photons enter the cavity, so that there is never more 
than one photon in the cavity at a time, and there is essentially no chance of a photon being 
coherently reabsorbed by the system from the detector. This last point is important for 
consistency. In that limit, we can trace out the atomic and reservoir degrees of freedom, 
and get an equation for the cavity mode alone, 



Pdot= -^detf^Pdet] + Tifepdet^ ~ y (^Pdet + Pd, 



ct 



+r 2 (b r b)p dct (tfb) - ^ ((Pb) 2 p dct + pdct(^) 2 ) • (3.3) 



2 

This equation is valid on timescales long compared to r. 
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In actual experiments the atoms are themselves measured as they leave the cavity, but 
this is unimportant for our present purpose. We can imagine them flying out into space 
forever, where they form (in the terminology of Gell-Mann and Hartle) a generalized record 
of the detector. Indeed, as we shall see in section V, the presence of the detector itself is 
not vital to derive a quantum trajectory description using decoherent histories. 

B. System and output mode 

Consider a quantum system with Hilbert space TCi, which is isolated except for a single 
channel of decay — an interaction with an external "output mode," which is continuously 
monitored. 7i 2 is the Hilbert space of the output mode — the cavity mode of ( |3.3|) . The 
combined state of the system plus output mode lies in the product Hilbert space 7ii <g> H 2 . 
This reduced model is illustrated in Figure 2. 

The measuring device produces two important effects. The first is dissipation. Excita- 
tions of the output mode will be absorbed by the measuring device at a rate Ti which we 
assume to be rapid compared to the interaction rate between the system and output mode. 
1/Ti limits the time-resolution of the detector. The second effect is decoherence. As the 
state of the output mode becomes correlated with the atom degrees of freedom in the detec- 
tor, the phase coherence between the ground and excited states of the output mode is lost. 
This loss of coherence is far quicker than the actual rate of energy loss. The decoherence 
rate is r 2 ^> IY 

Assume a linear interaction between the system and the output mode, and go to an 
interaction picture in the rotating wave approximation. This lets us remove the Hamiltonian 
of the output mode, and gives an interaction Hamiltonian 

Hj = K (a f ®b + a®ti). (3.4) 

The total Hamiltonian is 

H = H ® l + /t(a t ® b + a®b ] ), (3.5) 

where a and b (a^ and tf) are the lowering (raising) operators for 7i\ and 7i 2 , respectively. 
We want dissipation and decoherence to be rapid compared to the transfer rate between 
system and output mode. This requires that the system should not be too highly excited. 
The hierarchy of evolution rates is r 2 ^> T\ ^> K(a)a). We thus have a system with separated 
timescales. 

The system plus output mode obeys a Markovian master equation: 

p = -i[H, p] + r 1 SpS t - — {b% + pPb) 

+r 2 n 2 pn 2 - —{n 2 2 p + phi) = £p, (3.6) 

where p is the density matrix for the combined system and output mode, and £ is the 
Liouville superoperator. The operator h 2 = b^b acts on the output mode. Equation (|3.6|) is 
linear, and so can be formally solved: 
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p(tz)=exp(c(t 2 -t 1 )^p(t 1 ). (3.7) 

Since the output mode is heavily damped, we simplify the problem by retaining only the 
two lowest states |0) and |1), treating TC 2 as a two-level system. We then expand the density 
matrix p explicitly in terms of its components in Ji\ and H, 2 : 

pit) = poo(t) ® |0> <0| + poi(t) ® |0)(1| + pio(t) ® |1)(0| + pn(t) <g> (3.8) 

where the p^ are operators on Hi and the on Ti 2 . In terms of these components the 
master equation becomes 

Poo = -i[H , poo] - iKatpio + iKp 01 a + dpn, 

P01 = -i[H ,poi] - mcJpu + mpooa 1 " p i = p[ , 

pu = -z[# , P11] - ««apoi + mpioa 1 " - T 1 p n , (3.9) 

This model may seem highly simplified compared to an actual photodetection exper- 
iment, but it captures most of the essential physical principles without bogging down in 
unnecessary detail. To discuss photoemission, in is necessary to include some of the envi- 
ronment degrees of freedom explicitly. This is the function served by the output mode, which 
is about as simple as a physical system can be. The other important feature is the separation 
of timescales between the measuring device degrees of freedom and the effective dynamical 
timescale of the system. The usual approximation of ideal measurements as instantaneous 
depends on this separation of timescales. This sort of model has been used in the past to 
study the spectrum of emissions and time-correlation functions of an open optical system 



The important element in analyzing this model is its time evolution. Given that Ti <C T 2 , 
it is convenient to expand the time-evolution superoperator in the following form: 



e CSt = e C 2 8t + 



I dt'e^-^C^' 
Jo 

rSt rSt 

+ dt' dt"e^ 5t - t "^C 1 e c ^"- t ^C l e m ' + ■■■, (3.10) 
Jo Jf 



where multiplication of superoperators is composition, with the earliest rightmost. Second 
order terms are all that will be needed in this paper. Here the terms of the master equation 
have been separated: 

£ = £ 1 + £ 2 , (3.11) 

with 

dp = -i[H, p] + T 1 bpb j - —{b% + ptfb), (3.12) 



and 



C 2 p = T 2 n 2 pn 2 - —{n 2 2 p + phi). (3.13) 
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The effects of the superoperators £,\ and e C2t are given by 

(Ap)oo= -i[H , poo] - maVio + ^Poia + Tipn, 

(£ip)oi = -z[#o,Poi] - ina'pu + mpooa 1 " - yPoi = (Ap)io> 

(£ip)n = -«[#o> Pn] - ^apoi + mpioa 1 " - Iipii, 
(e £2i p)oo= Poo, 

(e £2i p)oi= (e- T2t/2 ) P oi = (e^p)L 

(e £2 V)n=Pn. (3.14) 

Since the superoperator e C2t is diagonal in the components of p, it is particularly simple 
to insert these expressions into the expansion (|3.10| ) and carry out the integrals explicitly. If 
we examine the evolution after a time 5t where Yi5t <C 1 <C Y 2 St, we see that the off-diagonal 
terms Poi,Pio are highly suppressed: 

2n 2 5t 

(e £5i p)oo~ Poo - i[H Q , poo}5t + Yip u 5t + — — (2a f piia - a f apoo - Pooa f a) 

1 2 

2 , . * . A + , 2i5t ~ . A . A+ . 2rx5t A . A+ . 
+ — (mpoia - wa'pio) - -p^[-Ho,««Poia - z^a'pioj — (map i - mpi a') 

1 2 1 2 J- 2 

-[#„, [#o, Poo]]^ 2 /2 - r lPll 5t 2 /2 - z[F , Y lPll }St 2 + O^/r 2 ), (3.15) 

2 , 2in8t * a , 

(e £<5 *p)oi ~ ^(^Pooa f -ma f pn) ^ — (-ia 1 "^, Pn] + «[#o, Poo]a f - ^i[a\pn]), (3.16) 

1 2 J- 2 



(e £<5t p)n~ Pn - i[#o,Pn]<ft - TiPiiSt + — — (2ap 00 a f - aa f pn - Piiaa f ) 

1 2 

2 A + . A 2i5t ~ . A + . A 2ri<5t,. A . A+ . 

+ — (z/tpioa 1 -mapoi) - -— |ii 0) mpioa 1 - map ij +— — (map i -iKp w a<) 

J- 2 1 2 1 2 

-[#„, [#o, Pn]]5t 2 /2 + r lPll 5t 2 /2 + i[^ , r lPll ]5t 2 + o^/r 2 ,). (3.17) 

The off-diagonal terms Poi, Pio will always be of order 0(k/Y 2 ). We can therefore consider 
an approximate set of differential equations in terms of poo and pn alone: 

Poo= -iH cS poo + ^Poo-^es + l^p u a + Y x p n , 

p u = -iHespu + ipuHls + japooa 1 - (Yi + j)pn, (3.18) 
where 7 = 4k 2 /Y 2 and the effective Hamiltonian 

H cS = H -i^a (3.19) 

is the same as that which appears in the quantum jumps formalism. This equation (]3.18|) 
is valid on timescales long compared to 1/T 2 . We examine it further in section HID below. 

If Yi is large compared to the other terms of the equation, then we can make the same 
sort of argument to show that pn will be highly suppressed (by a factor of roughly 7/T1) 
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compared to poo- in this limit we can therefore adiabatically eliminate all components other 
than poo HH- If we then consider a reduced density matrix p for the system alone, without 



the output mode, it will be essentially equal to p o (with a correction of order 7/T1 from 
P11). The equation for p then becomes 

p = —i[H , p] + ^apa* — ( r y/2)ci ! ap — ( , y/2)pa^a, (3.20) 

to first order in 7. This equation holds good on time scales long compared to I\. Thus, in 
the adiabatic limit we see that this indirect measurement scheme for the total system and 
output mode does reproduce the usual master equation ( |2.5|) for the system alone. Because 
7 is small, the damping is weak. This weakness is related to the quantum Zeno effect: if the 



detector had infinite time resolution, the system would never emit a photon at all E5] 



C. Quantum jumps and continuous measurements 

We can unravel equation (|3.20[) as described in section IIB into stochastic trajectories 



\ip(t,£)). By averaging over all realizations of with an appropriate probability 

measure (|2.10|) , one sees that it does reproduce the master equation (|3.20| ) as required ||. 



The master equation ( |3.20|) is valid only as long as the Markovian approximation remains 



good. In the case of our toy model, this means that it is valid only on timescales longer than 
1/Ti. Thus, rather than a jump occurring at a time ti, it is more correct to consider the jump 
as occurring during an interval At ~ 1 /Y\ centered on t^. For comparison to experiment this 
qualification is unimportant, but it will prove important in making comparisons to consistent 
histories. 

In the context of photon-counting experiments one can give a simple physical interpreta- 
tion to the individual quantum jump trajectories, as the state of the system conditioned on 
the continuous measurement record from the photon counter. As time passes without the 
detection of a photon we gain information about the state of the system; the lower states 
become more probable relative to the higher states, this effect given by the non-Hermitian 
part of the effective Hamiltonian. The jumps represent actual photon detections, in which 
both the state of the system and the state of our knowledge change abruptly. 

Consider the system plus output mode with Hamiltonian (|3.5|), and suppose that von 
Neumann measurements of the observable h 2 are performed repeatedly on the output mode, 
separated by time intervals 1 /T 2 . If we average over the two possible measurement outcomes, 
then the resulting mixed state is 

p^ApVa + ApA, (3.21) 

where Vq and V\ are projections onto the states |0), |1) of the output mode. The effect 
of this repeated measurement is to rapidly suppress the off-diagonal terms of the density 
operator p. 



The T 2 terms in the master equation (|3.6| ) have the same effect as ( 3.21 ) on time scales 



long compared to l/r 2 , both giving rise to (|3.18) in this limit. One can thus think of 



these terms being a continuous approximation to a series of repeated ideal measurements. 
However, terms of this form also arise generically in the study of systems interacting with 
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environments [Q. They are not unique to measurements. Indeed, this sort of description 
can be considered as a model of the measurement process itself: instead of taking place 
instantaneously, the measurement occurs over a short period of time, as the measured system 
becomes entangled with the many degrees of freedom of the measuring device. 

One must include as well the Ti terms which cause the emitted photons to be absorbed 
(either by the environment or by the measuring device, depending on the situation de- 
scribed). This could be modeled simply by having the repeated measurements reset the 
output mode to |0) after measuring it, with a rate of Ti (i.e., after every nth measure- 
ment where n = r 2 /T 1 ). The exact value of Ti is not very important. It represents the 
time-resolution of the measurement device, and all that is required is that photons not be 
emitted more rapidly than they can be absorbed. Indeed, if one merely assumes that the 
measurement device resets the state to |0) immediately after each measurement, but that 
the record cannot resolve individual clicks within an interval smaller than 1/Ti, the prob- 
ability of a given measurement record will be exactly given by ( 2.10| ); the evolution of the 
system conditioned on the measurement record follows the quantum jumps formalism given 
in section IIB. 



D. Ortho-jumps 



Consider the coupled equations ( [3.18D , which arise by averaging the full master equation 
( p.6| ) over times long compared to 1/T 2 . These equations are equivalent to a Lindblad master 
equation for the system plus output mode 

3 . . 1 . . 1 . . 

P = ~i[Ho, P]+J2 L ™P L m - 2 L ™ L mP - ^pLl.Lrn, (3.22) 
m=l 

with Lindblad operators 

Li = ^fya <g> ¥, 

L 3 = 1 <g> b. (3.23) 

Any density operator of the form p = poo|0)(0| + pn|l)(l| will remain of that form for all 
time. 

We can now unravel this master equation into the orthogonal jump trajectories of Diosi. 
For the moment, let us treat the deterministic and jump terms separately. Let us also 
continue to neglect the higher excitations of the output mode, so that it can be treated 
simply as a two-level system. If the system begins in an initial state = \<p) <8> |0), then 
the ortho-jump equation ( |2.19| ) becomes 

|#) = (-iH \4>) - |a f a|0) + ^(tfa)^)) <g) \0)dt + jumps. (3.24) 

If the system begins in the state = |0) <8> |1) then the equation is nearly identical: 

\dip) = (-iH \(f)) - |a f a|0) + ^(a + a)|0)J ® \l)dt + jumps. (3.25) 
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In both cases we see that the deterministic part of the evolution is the same on the system 
part of the state, leaving the output mode unchanged, and is equivalent to the deterministic 
part of the nonlinear quantum jump equation ( |2.13| ). What then are the effects of the jumps? 

A system in a state can jump into any of a set of orthogonal states which form a 
basis for the subspace spanned by the states \ipj) = (Lj — (Lj))\ip). If the system is initially 
in the state \<p) <S> |0) then only L\ makes a non-vanishing contribution to that space. Thus, 
there is only one type of jump from that initial condition, 

|0) <g> |0) -> a\<f>) ® \l)/(a ] a). (3.26) 

These jumps occur with a mean rate M\^(dNi) = ^(a)a)dt. Thus, the effect of this jump 
on the system is the same as a standard quantum jump, and they occur with the same rate. 

For a state |0) ® |1) the other two Lindblad operators L 2 and L 3 can give rise to jumps. 
However, because r x 3> 7, we can to an excellent approximation neglect the jumps due to 
L 2 , which occur very rarely compared to jumps due to L 3 . Making this approximation, the 
jumps take the form 

|$®|1)-|0>®|O), (3-27) 

and occur with a mean rate M(dN 3 ) = I\ independent of the state <fi. These downward 
jumps happen much more rapidly than the upward jumps due to Lj, and a trajectory has 
a near-unity probability of dropping back down to |0) within a time of order 1/Ti after 
jumping up to |1). 

We can thus characterize the typical evolution of these trajectories. They spend most 
of their time (by a ratio of roughly Ti/j) in the lower state |0), with the system evolving 
according to the deterministic terms of equation (2.13). At random intervals, there is a jump 



in which the system is multiplied by a and the output mode jumps to |1). Within a random 
time of order l/Ti the output mode drops back down to |0), while the system continues to 
evolve according to ( |2.13j) undisturbed. 

If we trace out the output mode, we see that one of these ortho-jump trajectories is 
exactly equivalent to a standard jump trajectory on the system Hilbert space alone. Each 
standard jump trajectory will have many ortho-jump trajectories equivalent to it, each corre- 
sponding to a slightly different time for the output mode to return to the state |0). However, 
as the exact time has no effect on the system evolution, this multiplicity is unimportant. The 
total probability of the ortho-jump trajectories will be exactly the same as the probability 
of the equivalent standard jump trajectory. 

This many-to-one relationship is a simple example of a coarse-graining; the ortho-jump 
trajectories are more detailed than the standard jump trajectories, because they include 
information about the state of the output mode as well. This concept is very important in 
the decoherent histories formulation of the same problem. What is more, it has already been 



shown by Paz and Zurek and Diosi |14| , |15| that orthogonal jump trajectories are equivalent 
to a particular set of decoherent histories. This new equivalence to standard jumps (true at 
least in this simple model) makes it easy to show how standard jumps can be represented 
by decoherent histories. 



Strictly speaking, the equivalence shown above in equations ( p.24| - pl25D is for the non 



linear version of both unravelings. As mentioned before, in general there is no linear version 
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of ortho-jumps. However, when neglecting jumps due to coherent reabsorption we see that 
this system is a special case, with all jumps between the orthogonal subspaces labeled by 
|0) and |1). 

Because of this there is a linear version of ortho-jumps for this case, produced simply by 
dropping all the nonlinear terms in equations ( |3.24| - pl25D and removing the renormalization 
in the jump (|3.26| ). The resulting equations are identical in form to the linear version of 
standard jumps, but the jumps are no longer strictly identified with "clicks" of a detector. 
The coarse-grained version of this linear equation is equivalent to standard linear quantum 
jumps on the system alone, just as in the nonlinear case. 



IV. CONSISTENT HISTORIES DESCRIPTION 



As described in section IID, a particular history is given by choosing one projection 
V^.(tj) at each time tj, specified by the sequence of indices {<x,} denoted h for short. The 
decoherence functional on a pair of histories h and h' is then given by 



D[h, h'\ = Tr Vl(t N ) ■ ■■V 1 ai (t 1 )p(to)V 1 a ,(t 1 ) ■ ~V2 N (t N ] 



(4.1) 



where p(to) is the initial density matrix of the system [|K| . 

We specialize to the system plus output mode described in section III. They are initially 
in the pure state |^) = \ipo) <S> |0). Since the degrees of freedom of the environment (e.g., the 
internal degrees of freedom of the measuring device) have already been traced out, we replace 



the simple Schrodinger evolution ( |2.26| ) with the Liouvillian evolution of master equation 
( |3. o] ), according to the quantum regression theorem [55]. The decoherence functional for 
two histories h and h! then has the form 



D[h,h'\ 



Ti{ V aN e 



est 



(V, 



e C5t (P ai \^)^\P ali )--.)P alN jP 



oin-V 



(4.2) 



where C is the Liouville superoperator from QSTEI). 

We consider histories composed of the following Schrodinger projections: 

p = i®|o)(o|, A = i®|i)(i|. (4.3) 



These projections represent the absence or presence of a photon in the output mode. These 
projections are spaced a short time 5t apart, and each history is composed of N projections, 
representing a total time T = N5t. A single history h is given by the string {a\, a2, . . . , ctjv}, 
where oij = 0, 1 represents whether or not a photon has been emitted at time tj = (j — l)St. 
The time-evolution superoperators in ( |4.2j ) tend to evolve pure states to mixed states. This 
is counteracted by the effect of the repeated projections P a , as we shall see. There are two 
important issues to address within the consistent histories formalism: the probabilities of 
histories (given by the diagonal terms of the decoherence functional) and the decoherence of 
the set of histories as a whole (given by the off-diagonal terms) . We examine them separately. 
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A. Probability of histories 



From the expressions (|3.10| - |37[4D , we can determine the character of the different histo- 



ries. The crucial choice is the size of the spacing St between projections. Too small and the 
histories will not decohere. Too large and all we will see will be standard master equation 
evolution, unresolved into trajectories. The interesting regime is in the range 

-!- « St « i (4.4) 

1 2 J- 1 



as described in equations ( p.l5| - j3TT8| ). On this timescale, the T 2 terms are sufficient to insure 



decoherence while the effects of the T 1 terms are resolved into individual trajectories. The 
time-evolution produced by the exp(CSt) superoperators is given by the simple equations 
( |3.18| ), equivalent to the averaged Lindblad equation ( |3.22j ). 



If the external mode is initially unexcited, with p = p 00 ® |0)(0|, then after evolving for 
a time St the state becomes 

( e CSt p\ 00 tt e -*(#o-i7a + a/2)5t^ ooe i(Ho+i7ata/2)<5t 



V)oi«o(«/r 2 ), 



a V)n~7aPooa f ^- (4.5) 



Here we see the appearance of the effective Hamiltonian H e a, just as in the quantum jump 
unraveling. 

We can also consider the case when the external mode is initially excited, with p = 
Pn <8> After a time St the state becomes 



(e £<5 V)oo~ TipuSt + ja^pnaSt, 
(e £<5 V)oi~ 0(k/T 2 ), 

(e £<5 V)n~ (1 - (ri + ^)St)e- l ^ St p n e l ^ St . (4.6) 

Once again the effective Hamiltonian appears, together with two additional effects. The first 
is the possibility that the photon in the excited mode will be absorbed by the measuring 
device. The second (much smaller) effect is the possibility that the photon will be coherently 
re-absorbed by the system. This process is so weak as to be negligible within the regime we 
are considering, and we will henceforth neglect it. This is the same approximation made in 
the ortho-jumps unraveling ( p.22| - [0^ ) by neglecting the L 2 contribution to the jumps. 



By combining the above expressions with the appropriate projections Vq and V\ (which 
pick out the poo or Pn component, respectively), we can write down the probabilities of 
different possible histories. Let us examine three illustrative cases and see how they exactly 
parallel quantum jump trajectories. 



1. Evolution without jumps 

Suppose that initially poo = \4>)(4>\ while poi = Pio = Pn = 0, i.e., the system is in a pure 
state and no photon has been emitted. Let us consider the history given by an unbroken 
string of iV P projections, corresponding to no photon being emitted during a time NSt. 
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The probability of such a history is given by the diagonal element D[0 N , N ] of (fO|). We 
can pick out the poo component of (4.5), and see that after the first time interval St 

V e cst (\ij)(ij\ ® |0)(0|)P « (e-^ ff<5t |^)(^|e^«' ata)5t ) ® |0)(0|. (4.7) 

Repeating this N times and tracing out the output mode Hilbert space 7i 2 we get 

£[0^,0^] « Tr! I e -^ eff ^ | ^> < V I e^e«^* 1 5 (4.8) 

which exactly agrees with the probability of the quantum jump or ortho-jump trajectories 
when no jumps are detected. (Tri indicates a trace over TCi only.) 



2. Evolution up to a single jump at time NSt 



Here we can make use of the previous result (|4.8| ) up until time NSt, when instead of 
using projection Vo we use V\. This is the same as keeping the pn component of exp(CSt)p 
instead of the poo component at the final projection time. This yields 

D[0"1,0*1] « (75t)Tr 1 |ae- i ^ ffAf<5t |V)(V'|e i ^ tffiV5 *a t }, (4.9) 

Once again, this agrees with the probability of the corresponding quantum jump tra- 
jectory. (Note, though, that in this context, St is short compared to the At in expression 
( p.lO| ). This is due to the finer-grained nature of this history, about which more below.) 



3. Evolution after a jump 



What happens after the external mode has "registered" as being in the excited state? 
Essentially, there are two possibilities: either the external mode can drop back down to the 
unexcited state (representing absorption of the photon by the measuring device) or it will 
remain in the excited state. We can examine these two possibilities separately: 

V e c5t (W)m (g) \1)(1\)A « T x 8tW){^\ ® |0)(0|, (4.10) 



Vie m (W)W\ ® |1)(1|)A « (1 - r 1 5t)e-^ ff<5t |V / )(^|e i ^ff <5 * g, ( 4 .11) 

So we see that the external mode has a probability of roughly T\St per time St of dropping 
back down to the ground state, whereupon it resumes evolution as in (|4.8| ), and a probability 
of 1 — TxSt of remaining in the excited state. In either case, the system component of the 
state continues to evolve according to the effective Hamiltonian H c q. 

This is exactly the same situation we encountered with ortho-jumps in section HID. 
There, we saw that a typical trajectory spent the majority of time with the output mode in 
the ground state. Infrequently this mode absorbs a photon from the system and becomes 
excited; this photon is dissipated into the environment in a time of order 1/Ti. The rate for 
coherent reabsorption by the system is much lower, and can be neglected. 
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Thus, we see that this set of histories corresponds exactly to the set of ortho-jump trajec- 
tories at the intermediate timescale St. This is, in fact, a special case of the correspondence 
shown by Diosi. In his treatment, the histories were branch- dependent in order to preserve 
the orthogonality of the jumps. In ours, the jumps are always between two known subspaces, 
and thus are branch independent. 

To reproduce standard jumps, we must coarse-grain just as before. Consider an interval 
At = M5t which is long compared to 1/Ti but still short compared to the dynamical 
timescales of the system. Let uj be a frequency which characterizes the system Hamiltonian 
H . Then we have TiAt ^> 1 ^> uj At ^> u5t. Let us sum all histories which include one jump 
and re-absorption within this interval (the probability of multiple jumps is small enough to 
be neglected), and look at the time-propagators occurring inside the trace: 

J2 Tri { ■ • ■ e - iH ^ jSt ae- lH ^ M -^ St — |VX^I — e lH ^ M - j)5t d^e lH ^ 5t • - • 1 

PS (M + l)Tn |. • . e -iH eS MSt/2 &e -iH eS M8t/2 . . . |^| . . . jit^ MSt/2 . . . j ? (412) 

with an error of order Ma; 2 At 2 . The exact absorption time is irrelevant, since the external 
mode is traced over, the system dynamics are essentially unaffected, and the probabilities 
sum to 1. So making this coarse-graining, and combining the above cases gives an expression 
for the probability of a particular jump record of 

p(ti, . . . , t n ) = ( 1 At) N Ti^e~ tAcs{T -' tN) ae' tAcs{tN ~ tN - l) a ■ ■ ■ ae~ iAcStl 

x|V)(V'|e^e* 1 a t • • • aV^e( T -*")\ (4.13) 

which exactly matches the quantum jump expression ( |2.10| ). 



B. Decoherence of histories 

Such a histories description is only meaningful if the histories decohere. Exact consis- 



tency, as in ( 2.22[) , is a difficult criterion to meet. It is more usual to show that a model is 
approximately consistent, which generally insures that the histories satisfy the probability 
sum rules to some level of precision. 

One criterion for approximate consistency has been suggested by Dowker and Halliwell 
||37|| . If we wish the probability sum rules to be satisfied to a precision e < 1, we require 
that 

\D[h,h'}\ 2 < e 2 D[h,h]D[h',h'} = e 2 p(h)p(ti), (4.14) 

for all unequal pairs of histories h, h'. Generally speaking, the "more different" a pair of 
histories is (i.e., the more projections they differ in), the more suppressed the off-diagonal 
term. This is certainly true for this model of photodetection. So it suffices to look at two 
histories which close as possible without being identical. 
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In the case of these "jump" histories, this means that these histories differ at a single time 
ti, one having a projection Vo, the other V\. In the decoherence functional, this is equivalent 
to picking out the p i or Pio component of exp(C5t)\ip') at that time, sandwiched between 
identical projectors V or V\ on either side. 

Examining the components given by ( |3.10| - |37L4] ) and ( |4.5| - |4.6| ), we see that 



\D[h,h'}\ 2 1 

1 1 (4.15) 



p(h)p(h') (T 2 5t) 2: 

so we expect the sum rules to be obeyed with a precision of roughly 0(l/T 2 St). For large 
T 2 this is more than adequate. 



C. Partial trace decoherence 

Finkelstein |38j has suggested an interesting alternative definition for decoherence in the 
case of a system-environment split, which he terms PT or "Partial Trace" decoherence. 
Consider the operator-valued functional defined by 

D[h, h'] = Tr eiDf \v^(t N ) • ■■Vl 1 {t 1 )p{t )Vl, i (t 1 ) ■ ■ -V^y (4.16) 

where Tr env denotes a partial trace over the environment degrees of freedom only. D[h, h'] 
is therefore an operator in the Hilbert space Tlx ® Ti.2- The criterion for PT decoherence is 

D[h,ti] =0,h^ti, (4.17) 

which is a much stronger condition than the usual decoherence criterion; it means that the 
different alternative histories are orthogonal in the environment degrees of freedom alone, 



implying the existence of generalized records [18 



From the form of the standard decoherence functional and the correspondence with the 
ortho-jumps unraveling (p.24^pT25|) , we see that D[h,h] will be a projector onto the state 



of the system plus output mode corresponding to the particular trajectory j. And from the 
equivalence of this to standard jumps, it follows that 

D[h,h] = $ h )p{h)$ h \ ®{\0)(0\ or |1)(1|), (4.18) 

where \iph) is the normalized state produced by the quantum jump trajectory corresponding 
to the history h, and p(h) is the probability of this trajectory. By exactly the same argument 
as before, the off-diagonal terms vanish and therefore this system is PT decoherent. 



V. GENERALIZED RECORDS 

As discussed in section IIB, one of the most powerful sources of decoherence is the 
creation of entanglement between "system" and "environment" degrees of freedom. Because 
these correlations exist, a careful measurement of the environment could in principle reveal 
the exact state of the system; thus, the different possible system states cannot interfere, 
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and condition ( |2.22| ) is assured. Gell-Mann and Hartle term these persistent correlations 
"generalized records;" they are generalized because in practice the measurement needed to 
access them might be impossibly difficult. 

For instance, suppose a measuring device records a series of results onto a paper tape. 
This is a record in the usual sense, and its existence guarantees decoherence between different 
possible measurement outcomes. If we then burn the tape, the record becomes difficult to 
access. But nevertheless, in principle the information is not lost; by cleverly measuring the 
ashes, air, and outgoing light one might in principle be able to reconstruct it, and thus 
decoherence persists. This "in principle" recoverable information is a generalized record. 

Suppose our system and environment begin in an initial pure state. Then the decoherence 
functional becomes (|2.25|) , implying that the unnormalized states {CV i |\I / )} are all orthogonal 
(or zero). This implies the existence of a set of orthogonal projection operators Rh such that 
RhCh'l^) = fihh'Chl^). Any such collection of orthogonal projection operators corresponds 
to an observable. Thus, in the pure state case, decoherence implies the existence of an 
observable whose value is exactly correlated with which history occurred. Conversely, the 
existence of such an observable implies decoherence. This converse holds even in the case 
of mixed states. This observable is the generalized record, and its existence is a stronger 
consistency criterion than ( |2.22| ). If the record is contained entirely in the environment 
degrees of freedom, so Rh = l sys ®-R^ nv where -R^ nv is a projector on the environment Hilbert 
space alone, it implies the partial trace decoherence (|4.17|) of Finkelstein. 

In the case of our model, the outgoing atoms clearly serve as such a generalized record. 
Suppose the output mode is in the ground state throughout some period At. In that case, 
all the outgoing atoms should be in the state (|0) + 1 1) ) / \/2- If the output mode is instead 
in the excited state, the probability is overwhelming that at least some of the atoms will be 
found in the state (|0) — |1))/a/2. Thus, these atoms serve as a record of the state of the 
output mode during the period At. 

Of course, if the output mode were only in the excited state for a short period, the atoms 
might fail to record the fact. But since photons in the output mode persist for an average 
time of order 1/Ti ^> At, this ambiguity is unimportant in practice. 

What is more, projections on the states of different atoms commute with each other. 
Thus we might, in principle, make no assertions about the state of the output mode, but 
merely wait until a large number of atoms had passed through the detector and then project 
onto their collective state. This projection would be sufficient to reconstruct the quantum 
trajectory in its entirety. 

In fact, we can play even more elaborate games than that. There is no reason, after all, 
why we should use the same projections on all the atoms. We might switch from one set to 
another as we chose. In the language of quantum trajectories, this would be like switching 
our choice of unraveling at different times. Indeed, we could go further than this: we could 
make our choice of projections on some atoms contingent on our results for other atoms. 
Thus, our choice of unraveling at one time might depend on the state at another time. What 
is more, there is no reason that we could not choose our unraveling at early times (i.e., the 
first atoms to pass through) contingent upon results at later times (i.e., the last atoms to 
pass through), rather than the other way around. Ordinary quantum trajectories theory 
does not consider this possibility for good reason, since such a measurement scheme would 
be extremely difficult in practice, nor is it obvious what is gained from it, but we can see 
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how the fairly specific assumptions of quantum trajectories can be progressively broadened 
towards the greater generality of consistent histories. The existence of generalized records 
insures that our set of histories will be consistent under all such permutations. 

The existence of generalized records also lets us get away from the whole notion of mea- 
surement. Suppose that instead of putting a photodetector outside our cavity we put nothing 
there, and allow emitted photons to escape to infinity. Away from the immediate vicinity 
of the system, we can describe this external electromagnetic fields in terms of incoming and 
outgoing plane waves with creation and annihilation operators a ±ki and a±kt, where £ labels 
the polarization. A reasonable initial condition would be that the incoming field is in the 
vacuum state. The outgoing field carries away information about the system. Using plane 
waves makes it difficult to discuss properties which are local in time; consider instead the 
operators c XQ k e and c\ QkQl defined by 

r ( (b — k ) 2 ) 

c Xo k e = (const) x J dkexp I -—^ h ikx > a M} (5.1) 

which annihilate (create) a Gaussian wave packet centered at xq and k . Needless to say, 
these states are overcomplete. But we can still use them to define operators Fc, which are 
quasiprojectors onto states with a single photon in phase space region C, and act as the 
identity outside region C. If C has a regular boundary, then as the volume of C becomes 
large compared to H Fc becomes closer and closer to a true projector. What is more, if C 



and C are nonoverlapping regions, Fc and Fqi approximately commute [^7 . 

Consider now a rectangular region C with length cAt much longer than the width of the 
wave packets, located well outside the near field of the system, and extent in wavenumber 
space Ik 3> Ak which is large compared to the linewidth of the emitted photons and limited 
to outgoing photons, centered on the cavity wavenumber ujq/c. Consider histories with 
projections Fc and 1 — Fc at a series of times tj = jAt. The projection 1—Fq is essentially a 
projection on no photon being emitted, since the probabilities of multiple photons or photons 
outside the frequency range are negligibly small. This set of histories will correspond on the 
system space to exactly the set of trajectories ( j4.13| ) we derived before. 



Time evolution by At is just equivalent to shifting the region C "outward" by a distance 
cAt. Thus the Heisenberg projections Fc{tj) which contribute to the history operator Q2.23|) 
are actually projectors onto non-overlapping regions, and hence commute. We could replace 
them all with the single projector 

K..i* = K(tn)---F h (t 1 ), (5.2) 

where ij = 0, 1 and Fi(tj) = F c {tj), F (tj) = 1 — F c (tj). The observable represented by the 
set of projections R is a generalized record of the quantum jumps trajectory. 

However, just as in the case of the photodetector, we can imagine projecting onto many 
other things besides photon number, each choice corresponding to a different unraveling. As 
long as we restrict our projections to outgoing waves, consistency is guaranteed. If we instead 
choose to project onto some combination of incoming and outgoing waves, the histories would 
not be consistent, in general. This is just as if our measurements of the outgoing field had 
a chance of reflecting light back into the system; in that case the condition ( |2.32 ) would 
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generally not hold, and the measurement scheme would no longer be an unraveling of the 
master equation. 

Again, we might choose one set of projections at one time and another at another, so 
as to change unravelings, and there is no reason in principle not to make the choice of 
projections at one time contingent on the results at another. 

This freedom leads to a somewhat curious situation. The outgoing field from the system 
constitutes a generalized record of the system evolution. But what is it a record of? By 
choosing different measurement schemes, or more generally different consistent descriptions, 
we can reconstruct very different-looking evolutions. It has long been known that on the 
one hand there were many different, incompatible unravelings of the master equation, and 
on the other many different incompatible sets of consistent histories. We now see that these 
can correspond to different, incompatible generalized records. 

VI. CONCLUSIONS 

We have seen how a continuous measurement can be described in terms of decoherent 
histories, and how the resulting histories match the corresponding quantum trajectory un- 
raveling of the master equation. The probabilities of the decoherent histories match the 
weights of the given trajectories, and the off-diagonal terms of the decoherence functional 
are highly suppressed, as we would expect from a measurement situation. 

Decoherent histories, as is widely known, must obey a condition ( |2.22| ) for their proba- 
bilities to have a consistent interpretation. What has been less widely appreciated is that 
quantum trajectories must obey a consistency condition (|2.32| ) of their own. The reason that 
this has received so little remark is that all commonly considered unravelings of the master 
equation automatically satisfy this requirement. These schemes rely on measurements of an 
outgoing field, which constitutes a generalized record of the system evolution, guaranteeing 
that both consistency conditions are obeyed. 

The existence of multiple incompatible records creates an interesting ambiguity in our 
description. Suppose that instead of measuring the field as it is emitted from the system, we 
wait for a long time and measure the field far away from the system. Our measurements then 
correspond to information about the system evolution long before, a sort of super-delayed- 
choice experiment. The choice of unraveling might be delayed until long after the evolution 
itself is complete. Were there no quantum trajectories until the measurement occurred? 
Does one trajectory "really" happen? Do all of them somehow coexist? 

These are the sort of interpretational questions upon which gallons of ink can be spent 
without putting an end to argument. From a practical point of view, it doesn't matter. The 
fact that these sets of histories are consistent means that we can choose an unraveling and 
make arguments based on the trajectories without fear that inconsistencies will invalidate 
our results. While there may be no actual measurement record, in the cases we consider the 
existence of generalized records serves the same function. While these generalized records 
might be hard to identify, and even harder to measure (especially if they depart at the speed 
of light), they insure that distinct alternatives remain distinct, and cannot subsequently 
interfere with each other. This is a valuable fact, especially since real measurements rarely 
approach the level of perfection commonly assumed for pure state unravelings. If these 
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unravelings decohere from each other, they can be treated like classical alternatives without 
fear of contradictions arising. Imperfect measurements can then be treated by further coarse- 
graining of these histories. 

What is more, if there are many different unravelings corresponding to different sets of 
decoherent histories, the choice of which to use can become a matter of convenience. In 
some situations a jump-like description is most useful; in others, a diffusion equation. It 
might be possible to use either one by an appropriate choice of projection operators on the 
same environment. 

Consistent histories is a powerful formalism, but producing a full set of histories can be 
dauntingly difficult, especially as their number increases exponentially with the number of 
projections. In practical calculations, it is often far easier to randomly choose some subset 
of "typical" histories and average over them. For this, quantum trajectory equations are 
supremely well-suited. Once again, the freedom to choose among different but equivalent 
formulations of a problem is a boon to the theorist. 

Over the last several years a great deal of work has been done on problems of quantum 
measurement, from a wide variety of viewpoints. I believe that these approaches, to the 
extent that they are correct, are connected at a deep level, and that by understanding the 
connections their power can be enhanced. This work is a small step in that direction. I have 
no doubt that much more progress is to come. 
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Figure 1. Picture of the atom-beam photodetector. The atoms are prepared in an 
initial state \tp) = (|0) + |l))/v / 2 and pass one at a time through the cavity where they 
interact briefly with the field of the cavity mode. The average time r between atoms is short 
compared to the lifetime of the cavity mode, which is in turn short compared to the rate at 
which photons enter the cavity. Thus, on average much less than one photon is present in 
the cavity. 



Figure 2. Schematic picture of the model system. The system, with internal Hamiltonian 
Ho interacts via an interaction Hamiltonian Hj with an output mode, which is modeled as a 
two-level system. The interaction is characterized by a rate k, and the output mode is part 
of a measuring device causing dissipation and decoherence with rates r!,r 2 . 
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Figure 1. Atom-beam photodetector. 
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Figure 2. Schematic setup of the model system. 



